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Recently, the decay properties of the double giant dipole resonance (DGDR) in several nuclei has been investigated 
experimentally in Coulomb excitation reactions at high energies [|^ . In particular, the neutron- and 7-decay channels 
were looked at in Ref. [^. In analysing the data, the authors rely on a model for the formation of the DGDR that 
involves the sequential excitation of the two phonon state through the one phonon state. Although it is concluded 
that the decay of the two phonons seem to follow the harmonic model (namely the two phonons decay independently 
from each other), we believe that a component essential to the analysis is missing, since the integrated cross sections 
obtained from the decay data deviate appreciably from the harmonic coupled channel calculations. The purpose of this 
paper is to develop a new model for the decay of the DGDR using the recently developed Direct -I- Fluctuation (DF) 
model of Ref. pM which reproduces the cross section value. We first give a brief description of the DF model. It is 
argued in Ref. M| that besides the, direct, g.s. —> one-phonon ^ two-phonon transition, there is another contribution 
that arises from the coupling of the one-phonon state to the complex background states followed by the excitation of 
a Brink-Axel phonon. The general structure of the two-step amplitude is 

Tfo = I Vf, (d + b) G. d V)o I (1) 

where Gi is the propagator in the region of the GDR, d and b are projectors for the one phonon collective states 
and for the background states responsible for the damping of the collective states respectively. In writing Eq. (|^) 
it is being assumed that the interaction I^o, when acting on the initial state, will not excite the background states 
appreciably. These states will however be reached through the damping mechanisms present in Gi, and will thus 
participate actively in the second step. We next think of the propagator Gi as split into a sum of average parts which 
are essentially diagonal in the subspaces b and d, Gb and Gd (with complex Q-value), plus an additional fluctuation 
part with zero average. Thus 
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and the average cross section becomes 
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with no cross term because of Eq. (^). 

The amplitude involving the averaged propagator is essentially the one given by the coupled channels calculation 
described in Ref. [^. The fluctuation part, Eq. (^, involves two different contributions. In the first the final state 
is reached from the intermediate one-phonon state, while in the second it is reached from the fine structure states to 
which the intermediate one-phonon state decays. The importance of this contribution stems from the possibility of 
collective dipole excitation of these fine structure states. On general grounds it may be expected to be proportional 
to the spreading width of the one “cold”-phonon states, and interference contributions with the two “cold”-phonon 
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amplitude to be small. This leads to an excitation cross-section to the energy region of the DGDR which is of the 
form 
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where the last term comes from T^q. This term is quite small compared to the first two and we drop it in what follows. 
It is known from work on giant dipole states in hot nuclei that the “hot”-phonon width T^ increases with excitation 
energy and then saturates. However, at E* = 2 E^^^ it is essentially equal to Moreover, for heavy nuclei 


pGGfli ^ pGBii^ Denoting ad^r = iTaoT 
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, we have for the integrated cross section 


CT' — crdir + < 7 /;. 
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When discussing neutron - or 7 - decay of the DGDR region, one has to multiply each of the two terms on the RHS 
of Eq. d) separately by the corresponding branching ratio. For neutron emission, experience has shown that besides 
the compound decay, there is also a preequilibrium component. Denoting the branching ratio by 

R-,f 

f 

, with T denoting the appropriate transmission coefficient for the n- decay of state i {d^ or d*) through compound or 
preequilibrium configuration (r), we can write for neutron emission. 

(Ed.) + B““P^ (Ed. )]+a [BP(t.n (Ed-) + (Ed-)]. (7) 

We next turn to the 7 - decay. It has been established that the 7 - decay of the GDR is composed of a direct 
plus compound pieces. Thus 


Bd^^ (Ed) = (Ed) + BH^P (Ed). (8) 

The decay of the DGDR proceeds in the following manner: the two dipole phonons decay directly in a sequential 
manner; B)]‘_i ,27 (Ed.) = , one phonon decay directly, while the other through the compound nuclear 

state at Ed; (Ed. ) = (Ed) (Ed), and finally the state simply mixes into the compound 

nucleus, which subsequently decays; (Ed.). Thus 

Bd2^27 (EdO = {Bi% (Ed))' + 2Bdi!^^ (Ed) B^^'’ (Ed) + B™ 2 T 27 ' (9) 


To simplify the form of Eq. ( 9 ), it in tempting to assume that (Ed.) = (Ed)^ , which would give 

for Bd2^27 (Ed.) = (Bd-,7 (Ed))^. However, for the moment, we will use Eq. ( 9 ) as it stands. 

The next step is to obtain an expression for the 7-decay of the d* state. It is a very simple matter to convince 
oneself that a collective Brink-Axel phonon first decays directly by emitting one 7, followed by the compound decay 
from the b space at Ed- There is also a compound-compound component, 


Bd.^27 (Ed.) = (Ed) Bn7'’ (Ed) + B“T27 (Ed.). 

Gollecting terms (Eqs. 9 and 10 ), we find for the 27 - emission cross section 
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+ crfi 


B^iU^ (Ed) Bn7" (Ed) + B™“P 7 (Ed. ) 


( 11 ) 


Equations ( 7 ) and ( 11 ) are the principal result of this work. In the following we present a detailed calculation for 
the DGDR decay in ^P®Pb. 

Galculations of the statistical decay of the giant resonance states were performed using the code STAPRE. This 
code permits the calculation of up to six sequential emissions from a compound nucleus, can calculate the complete 
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7 cascade for each residual nucleus in the decay chain and permits a simple pre-equilibrium decay calculation (no 
angular momentum or isospin dependence). 

The input data necessary for the calculations are the structural data of the residual nuclei, such as discrete state 
characteristics and level density parameters, as well as the transmission coefficients for particle emission and partial 
widths for gamma emission. In the calculations of the decay of described here, the level density parameters 

were estimated using a shell-corrected fit to their systematic dependence on mass and charge numbers. In particular, 
the level density parameters for ^°®Pb and ^°’^Pb were taken to be 0208 = 7-32 MeV“^ and 0207 = 8.48 MeV“^. The 
discrete levels used were taken from the ENSDF evaluated level files. [O We included the first 5 levels of ^°®Pb, up to 
E3,=3.71 MeV, and the first 7 levels of the neutron-emission residue Wb, up to Ea;=2.66 MeV. Neutron transmission 
coefficients were calculated using the global optical potential parameters of Becchetti and Greenless. [Q Charged 
particle emission was neglected due to the high Coulomb barrier and relatively low excitation energies involved. El 
gamma emission was calculated using the Brink-Axel approximation with the the ^°®Pb giant dipole resonance located 
at its observed position Ed=13.4 MeV with its observed width Td= 4: MeV. Higher multipole gamma emission was 
calculated using Weisskopf single-particle estimates. 

To study the effects of pre-equilibrium emission, the initial particle-hole configuration of each of the resonances 
was specified according to the standard model of their nature. Thus, to study the decay of the single giant dipole 
resonance, the initial excited nucleus population was taken to be in a 1 -particle-l-hole configuration with = 1 “ at 
an excitation energy of Ed=13.4 MeV. Similarly, the double giant dipole resonance was modeled as a 2-particle-2-hole 
configuration at E£;2=26.8 MeV, with 1/6 of the population having = O'*' and 5/6 having = 2+. The hot giant 
dipole was modeled as a 3-particle-3-hole configuration, in which one of the particle-hole pairs is attributable to the 
giant dipole resonance with the other two pairs being the lowest order contribution (in terms of particle-hole pairs) 
of the the ‘melted ’ giant dipole resonance which furnishs the hot background. The excitation energy E^;. and spin 
populations of the hot dipole resonance were taken to be the same as those of the double giant dipole resonance. 

The branching ratio obtained for single-step statistical decay of the single giant dipole resonance to the ground 
state of ^°®Pb is = 2.7 x 10“^. We found the branching ratio for occupation of the ^°®Pb ground state after 

the entire 7 cascade to be only about 10% larger, = 3.0 x 10“®. Thus, the single-step Brink-Axel El 

emission dominates the statistical 7 decay of the single giant dipole resonance. We note, however, that the direct 7 
decay branch is much more important than the statistical decay branch. The total branching ratio for 7 emission is 
actually several orders of magnitude greater than the statistical ratio we have calculated, Bd^^ = 1.9 x 10“^. 

We found the branching ratios for statistical decay of the double giant dipole resonance d^ and the hot dipole 
resonance d* to be very similar. In both cases, the single-step decay to the ground state is negligible, as would be 
expected given the initial spin populations (O'*" and 2+) and the dominance of the Brink-Axel El decay mechanism. 
The code STAPRE gives no direct information on double 7 emission, as it furnishes only the single-step occupations 
and the occupations after the complete cascade. However, we can infer a good deal of information from the cascade 
occupations. For both the double and the hot giant resonances, the branching ratio we obtained for occupation of 
the ^°®Pb ground state after the statistical cascade is ~ « 2.6 x 10“®. This is about 

3 times what we would estimate for the two-step emission process, which would give a branching ratio of about 

« 9.0 X 10“^°. However, a closer look at the occupation probabilities during the 7 cascade reveal 

that a good deal of the probability passes through the = 3“ state of ^°®Pb. This decay path necessarily involves a 
final E3 7 emission to reach the ground state and is not the path of interest to us. If we look instead at the occupation 
probability that decays directly from the continuum to the ground state, we get a better estimate of the double 7 
emission. For this (approximately 27 ) partial 7 cascade, we indeed find excellent agreement with our estimate based 
on the single giant resonance decay, 

B“ 2 T 2 , « B“T 2 , « (BH")^ « 1-0 X I0-® . (12) 

Using this relation, we can reduce the expression for the 27 -emission cross section, Eq. ( 0 ) , to 

(72^ iEd2) « ad^r (Bd^^ (Ed))^ + afi Bd-^^ (Ed) B™“p (Ed) (13) 

where we have made use of Eq. (|^). We now use the fact the the branching ratio for compound 7 emission is much 
smaller than that for direct emission, >> to neglect the fluctuation term, which yields, 

tT 27 (Ed^) W a dir (Bd ^7 (Ed))^ ■ (14) 

As the direct double excitation cross section adir is what has been called the harmonic term, we conclude that the 
27 -emission cross section should approximately agree with its harmonic estimate. We observe that this is consistent 
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with the experimental value of Ref. |^, j = 1.25(40), which within the experimental error, is equal to our 
value of 1. We emphasize, however, that the 27 -emission cross section is not proportional to the summed double and 
hot dipole resonance excitation cross section. As is clear from the discussion leading to Eq. the 27 -emission 

from the fluctuating term is greatly suppressed relative to that from the direct term. We thus expect the 27 -emission 
cross section to be proportional to the direct double excitation cross section, Odir-, alone. 

Turning now to the neutron emission, we first observe that the branching ratio for pre-equilibrium emission from 
the initial particle-hole configuration can be interpreted as the branching ratio for “direct” emission (as contrasted 
with statistical emission) from the giant resonance. Following this reasoning, we find that, of the 4.4% pre-equilibrium 
fraction of the neutron emission from the ^°®Pb single giant dipole resonance, 65% is emitted from the 1-particle-1-hole 
configuration. This yields a branching ratio of w 2.6 x 10“^ for “direct” neutron emission, which is of the 

same size as the branching ratio for direct 7 emission. However, it is extremely small compared to the branching ratio 
of the dominant statistical neutron emission. 

Pre-equilibrium neutron emission is substantially more important in the decay of the double and the hot giant 
resonances. We find that 20% of the neutron emission from the double giant dipole resonance is pre-equilibrium 
emission. For the hot giant dipole resonance, 15% of the neutron emission is pre-equilibrium emission. We can 
estimate the branching ratio for direct neutron emission from the double giant resonance by observing that, of the 
20% pre-equilibrium fraction of the neutron emission, 30% is emitted from the 2-particle-2-hole state. We thus obtain 
a branching ratio for “direct” neutron emission of w 6.0 x 10“^. We observe that this value of 6 % has the 

pleasing property of being close to the difference between the pre-equilibrium fractions for emission from the double 
and the hot giant dipole resonances. This lends force to its interpretation as the “direct” component of the neutron 
emission. For the case of the double giant dipole resonance, the remaining pre-equilibrium emission fraction of about 
15% corresponds to the emission that passes through the hot dipole resonance, that is, through the configurations in 
which at least one of the coherent dipole phonons has “melted” into the background. 

Independent of the division of neutron emission into pre-equilibrium and equilibrium components, the total branch¬ 
ing ratio for neutron emission is approximately the same for both the double and the hot giant dipole resonances, 
~ Bd*—>n- We can then write the neutron emission cross section, Eq. ( 0 ), as 

~ dir 4” ^//)Bd2_,n ■ (1^) 

The neutron emission cross section thus reflects the summed double and hot giant dipole excitation cross section. This 
result is also consistent with the experimental value of Ref. |^, = 1.06(35), which, within the experimental 

error, is equal to our estimate, <Jnl<yn°‘™ = {ddir + <^fi)/<^dir ~ 1.33, where the numerical value of the cross section 
ratio was taken from Ref. |^]. 

We conclude from our analysis that the 7 and neutron emission decay modes can make an observable distinction 
between the double and hot giant dipole resonances. In ^°®Pb, the 27 emission comes predominantly from the double 
dipole resonance, while neutron emission comes about equally from the double and hot dipole resonances. We also 
find that pre-equilibrium neutron emission from the double dipole resonance is slightly enhanced over that from the 
hot dipole resonance and, in both cases, is greatly enhanced when compared to the pre-equilibrium emission from the 
single giant dipole resonance. 


[1] See, e.g., H. Emling, Prog. Part. Nucl. Phys. 33 (1994) 729. 

[2] K. Boretzky et al., Phys. Lett. B 384 (1996) 30. 

[3] B.V. Carlson, L.F. Canto, S. Cruz-Barrios, M.S. Hussein, and A.F.R. de Toledo Piza, “Multiphonon and “hot”-phonon 
isovector electric-dipole excitations”, submitted for publication. 

[4] L.F. Canto, B.V. Carlson, S. Cruz-Barrios, M.S. Hussein, and A.F.R. de Toledo Piza, “Fluctuation contributions to double 
giant dipole resonance excitation cross sections”, to appear in Proceedings of the XX Brazilian Nuclear Physics Workshop, 
World Scientiflc, Sept. 1997. 

[5] D. Brink, Ph.D. Thesis, Oxford University, unpublished (1955); P. Axel, Phys. Rev. 126 (1962) 671. 

[6] C.A. Bertulani, L.F. Canto, M.S. Hussein and A.F.R. de Toledo Piza, Phys. Rev. C 53 (1996) 334. 

[7] R. Bonetti, M. Chadwick, P.F. Hodgson, B.V. Carlson and M.S. Hussein, Phys. Rep. 202 (1991) 171. 

[8] C.A.P. Ceneviva, N. Teruya, H. Dias and M.S. Hussein, Phys. Rev. C 55 (1997) 1246. 

[9] See, e.g., J.R. Beene et ah, Phys. Rev. C 41 (1990) 920. 

[10] M. Uhl and B. Strohmaier, STAPRE - A Computer Code for Partiele Induced Activation Cross Sections and Related 
Quantities, IRK Report IRK 76/01, 1976. 


4 



[11] J.K. Tuli, Evaluated Nuclear Structure Date File, Brookhaven National Laboratory Report BNL-NCS-51655, 1983. 

[12] F.D. Bechetti and G.W. Greenless, Phys. Rev. 182 (1969) 1196. 


5 



